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Abstract 

The evolution of the cosmological perturbations is studied in the context of 
the Randall- Sundrum brane world scenario, in which our universe is realized on a 
three-brane in the five dimensional Anti-de Sitter(AdS) spacetime. We develop a 
formalism to solve the coupled dynamics of the cosmological perturbations in the 
brane world and the gravitational wave in the AdS bulk. Using our formalism, 
the late time evolution of the cosmological scalar perturbations at any scales larger 
than the AdS curvature scale / is shown to be identical with the one obtained in the 
conventional 4D cosmology, provided the effect of heavy graviton modes may be 
neglected. Here the late time means the epoch when the Hubble horizon in the 
4D brane world is sufficiently larger than the AdS curvature scale /. If the inflation 
occurs sufficiently lower than the scalar temperature anisotropics in the Cosmic 
Microwave Background at large scales can be calculated using the constancy of the 
Bardeen parameter as is done in the 4D cosmology. The assumption of the result 
is that the effect of the massive graviton with mass me~"° > in the brane world 
is negligible, where e"" is the scale factor of the brane world. We also discuss the 
effect of these massive gravitons on the evolution of the perturbations. 
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1 Introduction and Summary 



Much attention has been paid for the possibihty that we are hving on a 3-brane in higher 
dimensional spacetime [|I], 0]. This brane world picture alters the conventional notion of the 
extra dimensions. Particularly if the bulk is Anti de Sitter (AdS) spacetime, the extra dimen- 
sions could be large or even infinite. The action describing the brane world picture is given 

by 

S = ^ J (i^Xv^^ (tZ^ + -Cr J d^Xy'-gbrane + j d!^ X^J - Q brane matter , (1-1) 

where TZ^ is the 5D Ricci scalar, I is the curvature radius of the AdS spacetime and = SvrG 
where G is the Newton constant in the 5D spacetime. The brane has tension a and the 
induced metric on the brane is denoted as gtrane- Matter is confined to the 4D brane world 
and is described by the Lagrangian Cmatter- We will assume Z2 symmetry across the brane. 

Recently, Randall and Sundrum (RS) constructed a simple model for a brane world 0. 
They assumed the effect of the matter confined to the brane is negligible compared with that 
of the surface tension. Their solution is described by the metric; 




rfs^ = - (dz^ - dr^ + 5ijdx'dx^) . (1.2) 



It has been shown that the usual 4D gravitational interactions are recovered on the 3-brane. 
One of the fascinating feature of their model is that the 5D spacetime is not necessarily 
compactified. 

In RS model, the 3-brane is Minkowski spacetime. Solutions for homogeneous expanding 
brane world are obtained by many people [4-14]. It has been shown that the evolution of the 
universe is identical with that of the conventional 4D cosmology at sufficiently low energies. 
However in the real world, the universe has inhomogeneity which leads to our structure of 
the universe [15-17]. This inhomogeneity can be observed today, for example, in the Cosmic 
Microwave Background Radiation (CMB). Then the cosmological perturbations in the brane 
world give direct tests for a viability of the brane world idea. In addition, the inhomogeneous 
fluctuations on the brane could be a powerful observable to probe the existence of the extra 
dimensions. This is because the inhomogeneous fluctuations on the brane inevitably produce 
the perturbations of the bulk geometry |T^. The perturbations in the bulk affect the motion 
of the brane in turn. Then, in general, the dynamics on the brane cannot be separated from 
the dynamics in the bulk. This could add a new property to the evolution of the cosmological 
perturbations and could reject the brane world idea. 
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To study the evolution of cosmological perturbations, we should treat the coupled system 
of brane-bulk dynamics. The problem has the similarity with the dynamics of the domain 



wall interacting with the gravitational wave, which has been investigated in 4D spacetime |T9 |. 
In our case, the matter on the brane is dynamical. This makes the problem very difficult. 
We should find a solution for the brane with the cosmological expansion and inhomogeneous 
fiuctuations. The most straightforward way is to solve the 5D Einstein equation, however, it 
would be difficult to carry out in general. 

In this paper we propose a new method to deal with the problem. We observe that the 
brane world cosmology can be constructed by cutting the perturbed AdS spacetime along the 
suitable slicing and gluing two copies of remaining spacetime. The point is as follows. If we 
choose a slicing to cut the 5D AdS spacetime, the jump of the extrinsic curvature along the 
slicing is determined. Because the jump of the extrinsic curvature should be equated with the 
matter localized on the brane, the matter on the brane is also determined. In other words, a 
solution for a brane with the given matter can be obtained by finding a suitable slicing. To 
find the suitable slicing for the given matter, we need two kinds of coordinate transformations. 
One is a large coordinate transformation which leads to the slicing which determines the 
background matter. Another is an infinitesimal coordinate transformation which leads to 
the slicing which determines matter perturbations. The coordinate transformations will be 
determined by imposing the conditions on the matter such as equation of state. More detailed 
procedures will be described in the next section. 

Our main result is 

^ = -2$o = const. (1.3) 
P 

at superhorizon scales and for late times when the Hubble scale H is sufficiently low {H ^ ^^^)- 
Here 6p is the density fiuctuations and $o is the metric perturbations in the longitudinal gauge 
in the brane world and we assumed the barotropic index of the matter is constant. The point 
to observe is that the solution ( |1.3D is identical with the one obtained in the 4D cosmology. We 
can also show that the late time evolution of the perturbations agrees with the one obtained 
in the 4D cosmology at subhorizon scales larger than the AdS curvature scale /. 

The assumption to obtain the above results is that the effect of the massive graviton with 
mass m > rrieff = /^^e"° in the brane world is negligible where is the scale factor of the 
brane world. We can understand the fact that the massive graviton with m > me/f can modify 
the evolution from the following arguments. For late times, the cosmological brane approaches 
to the RS brane. For the RS brane, the 4D gravity is recovered by the zero-mode of 5D 
graviton [3,20-23]. The Kaluza-Klein modes give the correction to the 4D gravity. However, 
in the Anti-de Sitter spacetime, the brane is protected from the Kaluza-Klein modes by the 
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potential barrier which arises from the curvature of the AdS spacetiem. For earher times, the 
cosmological brane is located at larger z in the coordinate (|1.2|) . The point is that, for larger 
2;, the potential barrier becomes lower. Then the relatively light graviton can interact with 
the brane. Thus for early times, the 4D cosmology will be susceptible to the Kaluza-Klein 
modes of 5D graviton. If the brane interacts with the 5D gravitational perturbations, the 
gravitational waves are inevitably emitted to the 5D bulk. It will cause the modification in the 
evolution of the perturbations. This picture is consistent with the result that the modes with 
large m > me// = /~^e"° can modify the evolution because vn^ff becomes smaller for earlier 
times. 

The paper is organized as follows. In section II, we describe our formalism in details and 
derive the background solution using it. In section III, we calculate the perturbations at 
superhorizon scales for late times H < l^^ using the formalism. In section VI, we calculate 
the perturbations at subhorizon scales for late times H < l^^. It will be shown that the 
evolution of the perturbations is identical with the one obtained in the 4D cosmology for any 
scales larger than the AdS curvature scale, if the effect of the massive graviton m > rrteff is 
negligible. In section V we study the effect of the massive graviton m > rrieff on the evolution 
of the perturbations. Finally we discuss the implication of our results on the brane world 
cosmology. In the Appendix, we listed useful formulas for calculations. 



2 Formalism 
2.1 Background 

We shall start with the Randall and Sundrum solution for the barne world 0. They considered 
a single brane with positive tension a in the 5D anti-de Sitter spacetime. Setting the surface 
tension of the brane by 

«:V=^, (2.1) 

and assuming the Z2 symmetry across the brane, they found a solution described by the metric 
( 11.21) . The brane is located at z = I (see Fig.l). From the metric (|1.2|) , we see the brane world 
is Minkowski spacetime. 

Next we will seek the brane world with the cosmological expansion. For this purpose, we 
note that the RS solution can be obtained by the following procedure. First cut the AdS 
spacetime along z = I and delete the AdS spacetime from z = I to the boundary z = 0. Next 
glue two copies of the remaining spacetime along z = I (see Fig.l). The jump of the extrinsic 
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Figure 1: Conformal diagram of the AdS spacetime. The thick hne shows the trajectory of the 
RS brane and the dotted hne shows the brane with cosmological expansion. The RS solution 
is obtained by deleting the AdS spacetime from z = I to the boundary z = (shaded region) 
and gluing two copies of the remaining spacetime. 

curvature at z = I should be equated with the matter at z = /. Thus we must put the suitable 
matter on the brane to glue the spacetime. Then the content of the matter on the brane is 
restricted as ( |2.1| ). The above argument implies that if we use the different slicing to cut the 
AdS spacetime, we need different matter to glue the spacetime. This is because the jump of 
the extrinsic curvature depends on the slicing we use to cut the AdS spacetime. Thus if we 
can find appropriate slicing to cut the AdS spacetime, we can put suitable matter resulting 
the cosmological expansion on the brane (see Fig.2). 

Now, we explicitly carry out the procedure to find the appropriate slicing. 
(l)Start with the AdS spacetime; 




(2.2) 
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Figure 2: The solution for a brane with cosmological expansion is obtained by deleting the AdS 
spacetime from y = to the boundary (shaded region) and gluing two copies of the remaining 
spacetime. 

(2)Make the coordinate transformation from the coordinate system {z,t,x^) to {y,t,x^) by 

z^l{f{u)-g{v)), r^l{f{u) + g{v)), (2.3) 

where u, v are the null coordinates of the new coordinate system; u = {t — y)/l,v = {t + y)/l 
and f{u) and g{v) are the arbitrary functions. The resulting metric is 



= - dt') + / , ,,, 5ijdx'dx\ (2.4) 

{fW-g{v)y {f{u)-g{v)y 



where f'{x) = df{x)/dx and g'{x) = dg{x)/dx. For the future convenience, we put 

M,t) ^ A f'i'^)9'{v) 2a{y,t) ^ ^ (2 Ci^ 

{f{u)-g{v)r {f{u)-g{v)r 



(3) Delete the AdS spacetime from y = to the boundary and glue two copies of remaining 
spcaetime. The jump of the extrinsic curvature at the brane (y = 0) is determined by the first 
derivative of the metric with respect to y; 

^ 1 fit /I) + g'{t/l) 
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.... ^ u rm+g'm i -f"mg'm + fmg"m \ 

' I [ fm - am 2 nt/i)g'{t/i) ) ' ^"-^^ 

where we denote the power series expansion near the brane as 

a{y,t) = ao{t) + ai{t)\y\ + (2.7) 

The jump of the extrinsic curvature should be equated with matter on the brane. Taking the 
5D energy momentum tensor as 

T^' = dtagi0,-p,p,p,p)6iy), (2.8) 
the junction condition can be read off as (see Appendix A) 



(4) Determine the matter content on the brane by imposing the equation of state 

p = wp. (2.10) 

Then it gives one constraint on / and g. There remains one freedom in / and g. Since 6^"*^^^°'*^ 
determines the time shcing in the brane world, it is a gauge freedom in the brane world. We 
fix the gauge degree of freedom by demanding that t is the cosmological time; 

-\fm-9my~ ^ ^ 



Combining ( ^loD and (|2ll| ), we can determine the function f{t) and g{t). The y dependence of 



the metric can be obtained automatically by replacing f{t) to f{u) and g{t) to g{v). Hence we 
obtained the coordinate transformation which leads to a brane with matter of given equation 
of state. 

The induced metric on the brane becomes 

ds'^ = -dt^ + e2°°(*^(5ijc/x^rfx^ (2.12) 



where 6"°'-*^ = {f{t/l) — g{t/l)) ^ is the scale factor of the brane world. From (|2.1CI| ) and ( p.llj ), 
we can verify that and p satisfy 



P + Mq{p + p) = 0, 

P H , 

3 36 ' 
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"0 = ^^P+ (2-13) 



where n^a = ASnG^. The former is the usual conservation of the energy. Since the term 
proportional to falls rapidly, the latter is identical with the Friedmann equation for late 
times. We show the solution of f{u) and g{v) for late times in Appendix A. The solution has 
two constants of the integration. We will normalize the scale factor as e""'^*"*'"'"'="*) = 1. 



2.2 Perturbations 

In the previous subsection, we obtained the brane world with the cosmological expansion. In 
the real world the universe has inhomogeneity which leads to our structure of the universe. 
Then it is important to obtain the solution for an inhomogeneous brane world. In this paper, 
we will concentrate our attention to the scalar perturbations. Unlike the homogeneous brane. 



we cannot place the inhomogeneous brane in the exact AdS spacetime as is shown in Ref |]18 
This is because the inhomogeneous perturbations in the brane world inevitably produces the 
perturbations in the geometry of the bulk. The perturbations in the bulk affects the motion 
of the brane in turn. Then the equations for metric perturbations and matter perturbations 
in the brane world can not be separated from the dynamics in the bulk. We should solve 
the 5D perturbations at the same time. The coupled equations for the brane dynamics and 
gravitational perturbations in the bulk are in general very difficult to deal with. 

The non-separable nature of the brane-bulk dynamics can be seen from the power series 
expansion of the 5D Einstein equation near the brane. We will denote the power series ex- 
pansion near the brane as in ( p.7|) . The dynamical variables of the brane are the potential 
perturbations ^q, curvature perturbations \l/o, density perturbations 6p and velocity pertur- 
bations V (see Appendix B). We have two equations from conservations of energy-momentum 
tensor T'^'^,, = 

((p + p)e°"t;)- = -3aoip + p)e''°v + 5p+ip + p)%, (2.14) 
and the trace part of the Einstein equation in the brane world from 5D Einstein equation; 

^0 + 4do^o + ao$o + 2(ao + 2d^)$o - ^e-2-o(2V2^o - V'^o) 



3 

In the conventional 4D cosmology, in addition to these equations, we have the equation 

$0-^0 = 0, (2.16) 
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for matter with no anisotropic stress. Then we have closed set of the equations. However in 
the brane world, the correspondent equation derived from the 5D Einstein equation is 

E2 = -e-2"«($o-^o + iVo), (2.17) 

where E is the non-diagonal comopnent and N is the {y,y) component of the metric 
perturbations. The equation contains E2, so we cannot have closed set of equations for the 
variables on the brane. This is because the inhomogeneous fluctuations on the brane inevitably 
produces the gravitational wave in the bulk, which gives the effective anistoropic stress to the 
perturbations. 

However, the procedure to obtain the background solution can be applied for the inhomo- 
geneous brane. The strategy is as follows. We first consider the perturbed 5D AdS spacetime 
in the coordinate system ( p.2|) . We assume the perturbations are small enough to treat by 
linear perturbations. Because there is no matter in the bulk, the perturbations should satisfy 
the vacuum wave equation in the bulk. In the coordinate system (|1.2|) , the wave equation 
can be solved easily with the help of the transverse-traceless (TT) gauge. In the 5D space- 
time, the free graviton has five independent components which include one scalar component. 
Thus there is one variable for the choice of the perturbed AdS spacetime. Next, we take the 
coordinate transformation (|2.3|) to provide the cosmological background. The transformation 
function f{u) and g{v) is determined by the background matter. The perturbations in the 
coordinate system {y, t, x^) is then easily obtained by the usual procedure of the coordinate 
t r ansf ormat ion . 

Once the perturbed AdS spacetime is obtained, one might attempt to cut the spacetime 
along y = and glue two copies of remaining spacetime as is done for the background spacetime. 
However, we need to be more careful. The presence of matter on the brane bends the brane. 
For the background matter we made coordinate transformation so that the brane is located 
at y = 0. However the matter perturbations also bends the brane. Then the perturbed brane 
is no longer located at y = (see Fig. 3). The perturbations evaluated at ?/ = is not the 
perturbations induced on the barne. Since the observers in the brane world are confined to 
the brane, we should evaluate the perturbations induced on the brane. Thus we should make 
(infinitesimal) coordinate transformation x^^ = x^ + ^'^^ to ensure that y = denotes the 
location of the brane. In general, the coordinate transformation makes Qyfiifi = t,x*) nonzero. 
These components can be gauged away, so we will take Qy^ = 0. Now the gluing can be 
performed so as to determine the cosmological perturbations. 

Then we impose the conditions on the matter perturbations to determine the variable for 
the choice of the perturbed AdS spacetime and the infinitesimal coordinate transformation. 
The point is that imposing the conditions on the matter is equivalent to solving the Einstein 
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Figure 3: Bending of the brane due to the matter perturbations. Due to the bending, the 
brane is not located at y = 0. Furthermore, the brane is not perpendicular to the y axis. We 
should make infinitesimal coordinate transformation. 



equation. Let us remember that we impose the equation of state p = wp ( |2.1CI|) on the matter 
in deriving the background solution, which gives the Friedmann equation and conservation of 
the energy on the brane (|2.13|). The same is true of the perturbations. 



We summarize the procedure to obtain the inhomogeneous brane. 



(1) Let us start with the AdS spacetime with linear scalar perturbations; 

rfs2 ^ fL\ Uz^ _ (1 + 20)dr2 + 2b,idx'dT + ((1 - 2^)5,j + 2^,^) dx'dx^) , (2.18) 



z 

where h^i denotes dhjdx^ and E^ij denotes d'^E jdx'^dxK The independent component of the 
scalar 5D graviton is one. We will use the transverse-traceless gauge to fix the gauge freedom. 
The gauge fixing conditions are given by 

0-3^ + V^^ = O, 
20 + V^6 = 0, 

6 + 2^-2V^^ = 0. (2.19) 

Using these conditions, the Einstein equation in the bulk becomes 

d'^h 3dh d^h 

TfT - -TT - TTT + = 0' 2.20 

OZ'' Z OZ OT'^ 
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where h = (f),b,'^ and E. Taking the solutions of the form h{z,T,x^) = {z/iyS{z)e *'^'^e 
the equation for S{z) is obtained as 

-Sizy+(m^-- 

z \ z 



S{z)" + ^S{zy +{m^- — ) S{z) = 0, =uj^- , (2.21) 



where S{zy denotes dS{z)/dz. The solutions are given by linear combinations of the Bessel 
function and the Neumann function J2{mz) +amN2{mz). The coefficient is determined by 
the the boundary conditions at z ^ oo. We take the boundary conditions that the positive 
frequency functions are ingoing at 2; — > cxo. Then the solution is given by 

h={jj^ j ^3 J dm hUk)HP{mz)e-'^^e^^^,Ch = ^,b,^,E), (2.22) 

where H2^^ is the Hunkel function of the first kind. From the gauge fixing conditions, the 
coefficients hm{k) satisfy 



2k^ + 3m2 
2k^ + 3m2 

E^{k) = PE^k), (2.23) 



bm{k) — 01 2 I o„ 2 Em{k)i 



where Em{k) is the arbitrary coefficient. This corresponds to the one degree of freedom of the 
5D scalar graviton and represents the spectrum of the gravitational waves emitted from the 
perturbed brane. 

(2) Make the coordinate transformation to provide the cosmological background; 

z = l{f{u)-g{v)), T = l{f{u) + g{v)), (2.24) 
The pertrubations in the cosmological background can be obtained using ( p. 51) as 

h = J e-2"(2^'*) I -0^ I dm /i„(fc)iJ^'V^e-"(^'*))e-'""(*'^)e'^^, (2.25) 

where J is the factor which comes from the Jacobian of the transformation (|2.24|) ; 

- = l{^-f\u)+g\v))=lae--, ^ = ^i^f^u) + g'{v)) = -la'c'^, 
cIt r)7 

— = lU\u)+g\v)) = la'e-, - = i^f ^u) - g\v)) = -lae^. (2.26) 
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In addition, we make (infinitesimal) coordinate transformations = + to ensure tliat 
tlie brane is located at y = 0. After imposing the gauge conditions Qy^ = 0, there remains one 
freedom of the coordinate transformation (^^). Then we take a slicing along the spacetime 
y = to cut the spacetime. 

(3) Cut the spacetime at y = and glue two copies of the remaining spcaetime along y = 0. 
From the junction conditions the matter on the brane is determined in terms of Em{k) and 

(4) Finally impose the two conditions on the matter perturbations and determine Em{k) and 
We will impose the condition on the anisotropic stress and equation of state of the matter. 

3 Evolution of Perturbations at Superhorizon Scales 

Following the formalism developed in the previous section, we calculate the evolution of the 
perturbations. To simplify the calculations, we first consider the long-wave perturbations in 
the brane world. We shall take the limit 



We will calculate the evolution for late times where the Hubble scales is sufficiently low 



We will take the assumption that the modes with me "° > Z ^ do not contribute to the 
perturbations in the brane world. Then we assume 



k 



0. 



(3.1) 



if </ 



-1 



(3.2) 



mZe-"o < 1. 



(3.3) 



The effect of these modes will be discussed in the section V. 



3.1 Calculation of Perturbations 



(1) In the k limit, we see only E survives in ( p.23| ). Then we shall start with 




(3.4) 
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(2) We make the (large) coordinate transformation ( |2.3| ). Since E does not change in this 
coordinate transformation, the metric is given by 

ds^ = e^P^y^'^ {dy^ - df) + e=^"(^'*) [6ij + 2^^^) dx'dx^. (3.5) 

Due to the bending of the brane by matter perturbations, the brane is not located at y = 0. 
We perform the (infinitesimal) coordinate transformation 

^x^+ e*", = (e, e*, e')- (3.6) 

After this coordinate transformation, the perturbed metric is given by 

ds^ = e^f^^y''^ (^{l + 2N)dy^ ~ {l + 2<l>)dt^ + 2Adtdy) 

^^2a{y,t) (^^(^^ _ 2^)^^^. ^ 2E,,j) dx'dx^ + 2B^idx'dt + 2G,idx'dy) , (3.7) 

where 

$ = ^^ + p'^y + p^\ 

^ = -a^-a'e, 

E = E + ^, 

B = i-e^^f-^'W 

A = e-e', 

N = e'+P'e + k'- (3.8) 

We will denote the power series expansion near the brane as 

$(y,t) = $o(t) + $i(t)|y| + ---. (3.9) 

We take the gauge condition G = A = and Bq = 0,Eo = (see appendix B-1). This 
determines ^* and ^ in terms of 

e = rdye + fo, fo = -e2-«^o, 
Jo 

^ = - r dye^^^-'^^e - Eo. (3.10) 
Jo 

Then we obtain the metric perturbations induced on the brane 

% = Pia+h 

No = ei+M, (3.11) 
13 



and the first derivative of tlie metric perturbations 

= e-2-o(-2a + 2«oa + 2(«i-/?i)^i)), 
El = Ei-e-2"«^^. (3.12) 



(3)We take the perturbed energy momentum in the 5D spacetime as 





6T 



M 
N 



/ 

\0 (p + p)e-"«t;,i SpSij 



(3.13) 



where we assume the anisotropic stress of the matter is zero. The jump of the first derivative 
of the metric perturbations should be equated with the matter perturbations on the brane. 
This junction condition is obtained as (see Appendix B) 



^1 

El 



-aiNQ + -K^5p, 





'5p ^ 6p^ 
-2{Pi - ai)e-"«i 



0. 



(3.14) 



Combining ( |3.12| ) and ( p.l4|) , we can write the perturbed energy momentum tensor in terms 
of e and E^^ik); 



K'^{p + p)e"°v 




-6 [aoio - "o^o + "1^0 
2 [^0 + 2«oa - (2ao + Sd^)^ + (A + 2«i)To 



El 



so- 



(3.15) 



(4)Imposing the constraints on the matter determines the unknown function C,^ and E^- First 
^0 is determined by the shareless condition; 

= e^"°A. (3.16) 
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The coefficient Em is determined demanding the equation of state Sp = clSp. In the coordinate 
(t, y), E can be written as 

E{t,y) = e-2"(s^'*) j dmH^^\mle-''^y^'^)l^Eme-'"'^^''^\ (3.17) 

Let us take the hmit m/e~" <^ 1. Then we can use the asymptotic form of the Hunkel function 
for small argument; 

HPiz)r^^ + ^ + Oiz^z'..), (3.18) 
where we neglect the over- all numerical coefficient. E can be evaluated as 

E{t,y) = jdm ( + le-2"(^'*)) /^^^^-.^.(i,,)^ (3^^g^ 

Using the Jacobian of the transformation 



dy 

dr 

dt 



= l{-nt/l)+g'{t/l))=laoe-''\ 

y=0 

= l{r{t/l)+g'{t/l)) = la,e-^\ 

S/=0 



dz 
dy 

dz 
'dt 



= -l{f\t/l)+g'{t/l)) = -la,e-''\ 

j/=0 



= l{f{t/l)-g'{t/l)) = -la,e-"\ 

y=0 



(3.20) 



we can verify the following equations about = e^"°^i and Tq = —c^^^Eq] 
«ia = Aa = fdm(-l + 2z^^)Eme-'^\ 

a^^y = [dm ( -2do - ime-°o + 2i "°^"° | E^e"™", 

a,^y = f dm (-4ao - 2dg + m^e"^"" - zmdoe-"° + 2z "° + + ^n^e"-"^", 



(3.21) 



and 



To = I dm [2- 2i^^ ) ^^e-*"^^ (3.22) 



2;2 72 _ , / r_rn2 



where we used do / , do/ ~ (HI) <^ 1. Then, we obtain 

n'^aiSp = J dm (6d^ + eimdoe""") ^^e"™^, 

K^ai^p = J dm (6wal - 6imdoe-°o + 2m'^e-^''°) ^mC"^""^, (3.23) 

where p = wp. We assume = w = cosnst. Then we observe that the equation of state 6p = 
w6p is satisfied if me~"° -C i/. It implies that Em should select the modes with me""" <^ H. 
Note that the assumption mle~°"^ -C 1 is consistent with the result that only the modes with 
me~"° < H contribute to the perturbations for late times H < l^^. 
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3.2 Evolution of Perturbations at Superhorizon Scales 



Let us evaluate the metric perturbations $o and \l/o in the brane world ( ^.11| ). From 
and ( |3:22|) , we obtain 

$0 = ^0 = J dmEme-'"'^ (3.24) 

For late times where the Hubble horizon is sufficiently larger than the curvature scale of the 
AdS spacetime H < r(t) is given by (see Appendix A) 

r{t) - /e-- (y) . (3.25) 

Because Em selects the modes with me~"" ^ H, we find that mr ^ 1. Then the metric 
fiuctuations are constant. The density fiuctuations becomes 

47rG45p = -3d^$o, (3.26) 

where we use k'^o" = ASirG^. We finally obtain the metric perturbations and density fiuctuations 
on the brane 

^ = -2$o = const. (3.27) 
P 

This is one of the main result of our work. We notice that the solutions for metric perturbations 
and density perturbations are identical with those obtained in the conventional 4D cosmology. 
Note that, if w c^, the modes with me~°'° > H should contribute to $o in order to satisfy 
6p = cl6p. Once these modes contribute to $0) "^o is no longer constant. This also agrees with 
the 4D cosmology. 



4 Evolution of Perturbations at Subhorizon Scales 

In this section we shall study the evolution of perturbations at subhorizon scales fce~"° > H. 
We derive the equation to determine Em{k) for each mode with k. We again use the assumption 
that the modes with me~"° > do not contribute to the perturbations in the brane world. 
Then we investigate the late time evolution of perturbations at subhorizon scales larger than 
the AdS curvature scale / > ke~"° > H). 
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4.1 Calculations of Perturbations 



(1) We start with the perturbed AdS spacetime 

ds'^ = Q (rfz^ _ ^ 2</))rfr2 + 2bidx'dT + ((1 - 2^)5ij + 2E^ij) dx'dx^) , (4.1) 
where 0, 6, ^ and ^ is given by ( |2.23| ). 

(2) The perturbations in (y, t, x*) coordinate is obtained by the coordinate transformation. 
The resulting metric becomes (see ( |2.25| ) and ( |2.26| )) 



ds^ = e^'^^^'"'^ [{1 + 2N)dy^ - {1 + 2^)dt^ + 2A dt dyj 

_^^2a{y,t) 2^-^^^. + 2E^ij) dx'dx^ + 2B,idx'dt + 2G^idx'dy) , (4.2) 



where 

$ = (/«')'e-'^0, 

N = -{lafe-^^(P, 
A = -2{faa')e-^^(P, 

G = {la)e-''b. (4.3) 
In addition, we perform the (infinitesimal) coordinate transformation by 

x^-x^+r, r = (e,e*,n. (4.4) 

We will take the gauge condition G = A = and Bq = 0,Eq = 0. This determines ^* and ^ in 
terms of as 



e = rdy{A + iy)+fo, fo = e2-o(5o-4 
Jo 

e = - rdy{G + e'^^-'^^e)-Eo. (4.5) 
Jo 

Then we obtain the metric perturbations on the brane 

iVo = No + ei+Pia, (4.6) 
17 



and the first derivative of tlie metric perturbations 

^1 = -tti^i - cxoio - «2^0 + ^1 - "0^0 - OilTo, 

B, = e-2"o(_24^ + 2doa-2(A-«i)To-io + e2"°5i-e2"o6o), 

El = Ei-e-'^^'^eo-Go. (4.7) 



(3) Combining tlie junction conditions ( |3.14| ) and (^.Tp , we can write matter in terms of ^'^ and 
K'dp = 2(a + 2doa + (2«2 + /?2-/??-2ai/3i)a 



-6 [oLQil + (a2 - ai/5i)^o - + "o^o + "1^0 - "lA^o 
2 (a + 2doa + (2«2 + 132- f3l - 2ai/3i)e 
+<l>i - 2^1 + io + 2«oio + (A + 2di)fo - (A + 2ai)iVo 



«:^(p + j9)e"°T; = 2^ - 2«oa - e'°°5i + 2(/3i - ai)To + e^"«Go + Aq, 







2e"^""^^ + 2^1 - 2Go- 



(4.8) 



(4)Imposing the condition on the anisotoropic stress ttt and the equation of state = dpjdp 
determines Q and Em{k). Substituting = e^^'^{Ei — Gq) into 5p — c^5p = 0, we obtain the 
following form of the equation for each k\ 

J dmF{t, m; k)Em{k) = 0. (4.9) 

Because the detailed form of F{t,m] k) is rather complicated, we omit it here. Defining the 
Fourier transformation of the function F{t, m; k) by 

F{t,m;k) = J dm'F{m',m;k)e-'"'\ (4.10) 

this equation becomes 

J dmF{m, m; k)E,n{k) = 0. (4.11) 

The problem is to find the eigenstate of the matrix F{m', m) with eigenvalue 0. Since F{t, m; k) 
is consisted from the combination of the oscillating function, the existence of the solution is 
very likely. 
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4.2 Evolution of Perturbations at Subhorizn Scales 



The equation to determine Em{k) is rather comphcated. However we can deduce the evolution 
of perturbations for late times H < from the following arguments. We take the assumption 
that the massive modes with me""" > l^^ can be neglected. We can evaluate the metric 
perturbations as 



2k' + 3m2 

3m2 1 



lUJT 

1 



where we neglect the terms of the order (9((m/e""°)^). Then we find that 

$0 - ^0 = C((fc/e-"»)2). (4.13) 

For kle~"° but not necessarily fce~"° < H, these agree with the result obtained in ( p.24| ). 
A notable point is that this equation is sufficient to show that the 4D cosmology is reproduced. 
From the 5D Einstein equation, we have alredy obtained the three equations for $0) 
and V (conservation of the energy momentum ( |2.14| ) and trace-part of the Einstein equation 
( p.l5| )). Thus, for kle'""" -C 1, we have closed set of equations about the metric perturbations 
and density fluctuations, which is identical with the one obtained in the conventional 4D 
equation. Then we conclude that if the effect of the massive modes with m/e""" > 1 can be 
neglected, the cosmological perturbations are not interfered by the extra dimension. Whether 
the modes with m/e^"" > 1 contribute to the perturbations or not is determined by Em{k) 



which can be obtained from (|4.11|). 



5 Deviation from 4D cosmology. 

In the previous two sections, we show that the late time evolution of the perturbations is not 
modified by the 5D graviton if the the effect of the modes with me~"° 3> is negligible. 
Let us investigate the effect of these massive modes. To simplify the calculation, we assume 
k 0. Without taking the limit mle~'^° <^ 1, ai^o ^"^^ o^o^ becomes 

ai^^ = - J dm{mle-"")E^Hi^\mle-"°)e-''^^ 

+doe""° J dm{imf)EmH!^\mle-''°)e-''^\ 

aofo = alP J dm{mle-'''')E^H'i^\mle-'"')e-'"'^ 

-doe-"« J dm{imf)EraH^^\mle~'^'')e-''^\ (5.1) 
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Then we can evaluate \l/o as 

^0 = (1 - «o'/') / dm^Ut)Em{k)e-'''^, ^^(t) = (m/e-"°)i/f )(m/e-"°). (5.2) 

We shall investigate the effect of the massive modes. The evolution of can be estimated 
using the asymptotic form of the Hunkel function for large argument h[^\z) oc Xj \fz. We find 
that ^to(^) behaves as 

oc e-"«W/^ (m>/"^e"«). (5.3) 

Then, the modes with 

m > me// = C/p = (5.4) 

modify the evolution of the metric perturbations and hence of the density fluctuations. Because 
we have normalized the scale factor as e"o(*=*2"'='"="*) = 1, Ze// becomes larger than / for early 
times. 

The result can be interpreted as follows. For late times, the trajectory of the brane is 
identical with the RS brane (see Fig. 2). Thus the behavior of the gravity in the brane world 
for late times can be deduced from the arguments for RS solution. In the coordinate system 



(2;,r, x*), the perturbations h{z,T,x^) satisfies the wave equation (|2.2CI|) . Defining h{z,T,x'^) = 

wave equation can be rewritten as 

(-^^ + W^))^W = ^^¥(^), K4.5(^) = y^. (5.5) 

The problem can be understood as the potential problem in one dimensional quantum me- 
chanics where m represents the energy. If we cut the AdS spacetime aX z = I and put the brane 
there, the potential term proportional to 5{z — I) appears. Imposing the Z2 symmetry across 
the brane, the potential becomes a "volcano" potential. The term proportional to the delta 
function is responsible for the normalizable zero-mode of the 5D graviton which reproduces the 
4D gravity on the brane. The brane is protected from the the massive modes by the potential 
barrier Vacis (see Fig. 4). Only the modes with large m can affect the gravity on the brane. 
Let us consider the cosmological brane for early times. For early times, the brane is located 
at large z > I. For larger z, the potential barrier VAds is lower, then the mode with smaller m 
can affect the gravity on the brane. This picture is consistent with the result that the modes 
with large m > rrieff = /"^e"" can modify the evolution of the metric perturbations because 
rrieff becomes smaller for early times. 



6 Discussions 



The cosmological perturbations in the brane world provide useful tests for the brane world 
idea. This is because the perturbations in the brane world interact with the perturbations in 
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Figure 4: Volcano potential for 5D gravity. The location of the brane determines the height 
of the "barrier" which protects the brane from massive modes. 

the bulk which is inherent nature of perturbations in the brane world. The dynamics of the 
brane can not be separated from the dynamics of the bulk. This is because the inhomogeneous 
fluctuations on the brane inevitably produces the gravitational waves in the bulk, which in 
turn affect the evolution of perturbations on the brane. Thus, naively, we think the evolution 
of the cosmological perturbations is modified significantly. 

We showed that this is not a case for late time H < l^^ and at scales larger than the 
AdS curvature scale. The metric perturbations become frozen once the perturbations exit the 
horizon as in the conventionoal 4D cosmology. This result is important for the inflationary 
scenario in the brane world. If the scale is sufficiently higher than the scales of the inflation, 
and if heavy graviton modes may be neglected the constancy of the curvature perturbations 
can be used to estimate the scalar temperature anisotropics of the CMB at large scales. Our 
results are consistent with those of Ref [^, where curvature perturbations on large scales is 



shown to be conserved, and the density perturbations generated during high-energy inflation 
on the brane are calculated. 

The assumption in obtaining the above results is that the effect of massive gravitons with 
m > rrieff = e"°l~^ can be neglected. The contribution of these modes depends on the initial 
spectrum of the fluctuations. If the primordial fluctuations are generated during inflation at 
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low energies H < I ^, heavy gravitions are significantly suppressed p6[. Then, in this case, 
the assumption seems to be natural. 

Another modification of the evolution arises if the scales of the perturbations becomes 
smaller than /. This modification is not important for late time evolution because the cosmo- 
logical scales is significantly larger than /. However at the beginning of the infiation, we should 
consider the scales comparable to /. Then the modification becomes important to predict the 
primordial spectrum of the fiuctuations during inflation. 

The key to quantify these modifications is the understanding of Em{k). The equation for 
Ejn{k) obtained in this paper might give a way to estimate the effect of the interaction with 
bulk graviton on the perturbations in the brane world, which is the intrinsic feature of the 
brane world cosmology. 

Finally, we comment on the possibility for the extension of the present work. In this paper, 
we have considered the spatially flat universe. The extension to the non-flat universe may be 



possible using the coordinate system in Ref WM. The extension of our method to vector and 



tensor perturbations is straightforward. The tensor perturbations were calculated in the de 
Sitter background and the agreement with 4D gravity was demonstrated in a certain limit . 
It will be interesting to investigate the effect of the vector components of 5D graviton on the 
4D cosmological perturbations. 



While the present work was being completed, related papers Ref [26-29] appear on the 
hep-th. 
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A Background - Einstein equation and junction condi- 
tion 



In this appendix, we derive the junction conditions and obtain a solution for the background. 
The junction conditions can be obtained from the 5D Einstein equation; 

G'^ = ^,6^' + K'^^^^n^ = p'^ + e-^K'n^, iM,N = y,t,x^). (A.l) 

We take for the energy momentum tensor in the 5D spacetime 

= {-adtag{0, 1, 1, 1, 1) + diag{0, -p,p,p,p)) 5{y). (A.2) 

The jump of the first derivative of a{y,t) and P{y,t) gives the 6{y) function to the Einstein 
tensor. The Einstein tensor is given by 

G° = -3e-2^(/?a + a'/5-d'-da'), 

G) = 5]e-^'^{-2a-2d? -'l3 + 2a" + 2a'^ + (3"). (A.3) 
Equating the coefficient of 5{y) we obtain the junction conditions 

Following Ref we make power series expansion of the Einstein equation near the brane. 
The y^ order of the (y, y) ,{y, 0) component of the Einstein equation gives 



ao + 2al = 47rG4 (^^ - 



K,'^p{p + 3p) 



36 

p + 3do(p + p) =0, (A.5) 
where K^a = 487rG4. The integration of the first equation gives 

al=^-^p+^ + e-'«C, (A.6) 

where C is the constant of the integration and is proportional to the mass of the AdS- 
Schwarzshild mass. Because the bulk is AdS spaceatime in our solution, G = 0. These 
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equations are equivalent with ( p.l3|) . The order of (0,0) and components of the 



Einstein equation give ^2 and (32 in terms of ao, ai, (3q and 

2 

Ci2 = + "o/5o - + ai/^i + — , 

2 

/32 = d^ + 2ao + /5o + ai-2ai/5i + -. (A.7) 

Let us find the function f{u) and g{v) for late times. For late times, we can neglect 
0{{p/ (jY). Then the solution for e~"° is given by 

2 

e-"o = /(t//) - git /I) = a, (y)^™ ■ (A.8) 
Combining this equation with e^'^" = 1; 

, f'mg'm 



U{t/l)-g{t/l)Y ' 
we can obtain the first order differential equation for f{t/l); 



t\ 3(1+™) / / /S ^ t 



(A.9) 



For late times t/l ^ 1, we obtain the solution for f{t) and g(t); 

where a,, and 6* is the constants of integration. Then we obtain 

r(t) ~ /e-"» , (A.12) 

where we take 6* = 0. 



B Perturbations - Einstein equation and junction con- 
ditions 

B.l Gauge fixing and Einstein tensor 

In this section, we derive the 5D Einstein equation for perturbed AdS spacetime and derive the 
junctions conditions. We will concentrate on the scalar perturbations. We put the perturbed 
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5D AdS spacetime as 

ds^ = e2^(^'*) [{1 + 2N)dy^ - {1 + 2<i>)df + 2 A dt dy) 

+ e^"^^'*) (((1 - 2^)5ij + 2E^ij) dx'dx^ + 2B^idx'dt + 2G,idx'dy) . (B.l) 

There are three degrees of freedom in the gauge transformations 

By this gauge transformations, metric perturbations are transformed as 





= + + + 






E 


= E + C, 


B 




A 




G 




N 





Using these degree of the freedom, we impose the gauge conditions so that the resuhing 
coordinate becomes Gaussian Normal (GN) coordinate because the metric perturbations in 
the GN coordinate on the brane is the metric perturbations observed by the observers confined 
to the brane. In the GN coordinate, the transverse component of the metric gy^, {/j, = t, x*) 
vanish {G{y,t,x'^) = A(y,t,x^) = 0) and the brane is located at y = 0. The former conditions 
are achieved by ^ and ^* and the latter condition is achieved by The conditions G — A — 
determines ^* and ^ as 

e = l^'dy{A + e)+e\t,x'), 

e = -£dy{G + e'^f'-'^^e)+e{t,x'), (B.4) 

where e* and e are functions with no y dependence. These residual gauge transformations 
enable us to impose two additional gauge fixing conditions. We take -Bo(y = 0, t, x^) — EQ{y — 
0, t, x') = gauge on the analogy of the longitudinal gauge in the conventional 4D cosmological 
perturbations theory. The Einstein tensor is calculated as 

SG'^o = e"^^ ((6d^ + 6d/^)$ - 3a'N' - 3aN - {6a" + 12a^ - 6a'p')N 
-3^" + 30 + 2d)^ + 3{P' - 4a')*') 

{e-^f^{E" -($ + 2a)E - ((3' - Aa')E') - 2e-2"* + e'^^Ar 
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5G\ = e~^^(-(6a'^ + 6a/30A^ + 3a'$' + 3Q;6 + (6a + 12Q;^-6Q;/3)$ 

+3^ - 3(/5 - 4a)^ - 3(/?' + 2a')^') 

(e-2^(-E + (/3 - 4d)^ + + 2a')E') - 26-^"^ + e'^"* 

+e-2(a+/3)(^^(2d-/3))s), 
^G^o = e-2'^(3^' + 3d$' + 3a'iV-3(/5'-a')^-3(/9-a)^') 

(JG^ = e-'^ ((2q;' + /3')A^ + (a' - - ^' + 2*' 

+e-^^ (^(3da + a - 2a'/3)5 + (/? - ^a)B' + a'E - ^5 , 

5G^. = e-^^(-(2d; + /3)$-(d-/3)Ar + 7V-2* 

+e-2^ (-(3a'2 + a" - 2a'/3')^ - - ^a')^' + \b")) , 

5G) = e-2'^(5^($''+(4d + 2/3 + 6d2)$ + (/j + 2Q;)6+(/3' + 2Q;')*'-A^ 

-(Aa" + 2/3" + 6q;'^)A^ - + 2q;')A^' - (/3 + 2a)N + 2* + Gd* - 6q;'*' - 2*") 

(^e"^"($ + - ^) + e-^f^ (^e-^"{B + d^) - 3aE + 3a' E' -E + E"^^ 6) 
- (^e-^"($ + N-^) + e-^^ (|e-^°^(S + dS) - Sd^; + 3a' E' -E + ^")) * , (B.5) 
where B — e^"S and we denote h' — dh/dy and h — dh/dt. 

B.2 Einstein equation and junction conditions 

In this subsection, we derive the junction conditions. The Einstein equation is given by 

5G^ = K^e-^iSTHf - NT^). (B.6) 

We take for the 5D energy-momentum tensor 

/ 

-5p -{p + p)e^-v, I 5{y). (B.7) 
V (p + p)e-"owi 5p5ij 



jrrpM 



By equating the coefficients of S{y) in the Einstein equation gives the following junction con- 
ditions 

6 
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E, = 0, (B.8) 



where we use = 1. 
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